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UserGroup 
1993 U.K. Meeting 

in the frame of 

The Technology in Mathematics Conference 

at 

The University of Birmingham 

Monday, 20th September 1993 
10.00 a.m. till 4.00 p.m. 

**************************************************** 

ATM ATM ATM ATM ATM ATM ATM ATM 

In the autumn this year MicroMath will be running a special feature on the use of computer algebra 
systems and we have a number of articles specifically about Derive. Non- ATM members who would 
like to receive a copy should contact Basil Blackwell Publishers, 108 Cowley Road, Oxford, England. 

ATM ATM ATM ATM ATM ATM ATM ATM 

Announcing a NEW journal: 

The International DERIVE Journal, 

Reporting on Research and Innovative Teaching 

The International DERIVE Journal is published to disseminate information about research and practice 
on the use of DERIVE as a tool for doing and learning mathematics. The journal aims to enhance the 
use of DERIVE by reporting on research and significant innovations. 

Editor: Prof. John Berry, Centre for Teaching Mathematics, Plymouth, UK 


DERIVE - BOOK - SHELF 


Two Italian DERIVE books: 

Manara-Perotti, Lineare a Geometria con DERIVE, McGraw-Hill, Milano 1992 

Bacchelli-Lorenzi-Perotti, Matematica con DERIVE, McGraw-Hill, Milano 1992 

Mauve/Moos, Mathematik mit DERIVE (Arbeitsblatter zur experimentellen Mathematik), 

Heidelberg 1993 

***************************************************************** 
Ein Service der DUG fur den deutschsprachigen Raum: 

Alle DUG-Mitglieder konnen den Tagungsband der DERIVE-Konferenz 1992 
"Teaching Mathematics with DERIVE" zum Preis von oS 300.- (Oster- 
reich) , bzw. DM 45.- (Deutschland, Schweiz, ....) beim Herausgeber 
beziehen. Ein Anruf Oder eine kurze Nachricht geniigt. 
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The Derive-News-Letter is the 
Bulletin of the Derive-User 
Group. It is published at least 
three times a year with a con- 
tents of 30 pages minimum. The 
goals of the D-N-L are to en- 
able the exchange of experi- 
ences made with Derive as well 
as to create a group to discuss 
the possibilities of new me- 
thodical and didactic manners 
in teaching Mathematics. 

Editor : Mag. Josef Bohm 
A-3042 Wurmla 
D ' Lust 1 
Austria 

Phone/FAX: 43- ( 0 ) 2275/8207 
e-mail : no j o . boehmOpgv . at 


Contributions : 

Please send all contributions 
to the above address. Non- 
English speakers are encouraged 
to write their contributions in 
English to reinforce the inter- 
national touch of the D-N-L. It 
must be said, though, that non- 
English articles will be warmly 
welcomed nonetheless. Your con- 
tributions will be edited but 
not assessed. By submitting ar- 
ticles the author gives his 
consent for reprinting it in 
the D-N-L. The more contribu- 
tions you will send to the Edi- 
tor, the more lively and richer 
in contents the Derive-News- 
Letter will be. 


Preview: (Contributions for the next issues) 

Bisection with DERIVE, D.M.Dyer, USA 

Fluid Flow in DERIVE, Reuther a.o., BRA 

Newton's Chaos, Graphic Integration, J. Bohm, AUT 

Computer Aided Mathematics in School, K.H.Keunecke, GER 

Continued Fractions, R. Setif, FRA 

DERIVE in Hawaiian Classrooms, E. Sawada, USA 

Stability of Systems of ODEs, A. Kozubik, Slovakia 

Minimization of a „Flat Function 44 , Lopez, Portugal 

Digital Filter Design Using DERIVE, Hood, UK 


(DNL#1 1 will be published September 1993) 

Impressum : 

Medieninhaber : DERIVE User Group, A-3042 Wurmla, D'Lust 1, AUSTRIA 
Richtung: Fachzeitschrif t 
Herausgeber: Mag. Josef Bohm 
Herstellung: Selbstverlag 
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M. Hornschuh, Hilter, D 

Ich suche DERIVE - User im Raum Osnabriick zum Erfahrungsaustausch. 
Meine Telefonnummer: 05424 69370 


P. Steiner, Klagenfurt, A 

Can you tell me, how DERIVE calculates the values for SIN(x), COS(x), EXP(x) and LN(x)? 

DNL: I hope to get an answer from Hawaii. 

R. Mauve, Hirschhorn, D 

1 . Komplexere Grafiken im R 2 werden erst geplottet, wenn man die, sie erzeugenden Listen verein- 
facht hat (Simplify oder approX). Das fuhrt zu erheblichen Problemen bei der Erzeugung umfang- 
reicherer Zeichnungen. Z.B. zeichnet bei geeigneter Definition von VERSCHIEBE, RECHTECK 

VECTOR (VERSCHIEBE (RECHTECK ([1,1], [4,6]), [n,0]),n,-6,6,2) 

trotz einiger "Bemiihungen" nichts. Aber nach Vereinfachung mit Simplify ergibt sich: 



Wunsch: DERIVE sollte ohne Vereinfachen ("unausgerechnete Ausdriicke") plotten. 

(Wish: DERIVE should plot directly without Simplifying the expressions; it is not able to 
plot complex graphs without simplifying or approximating the generating lists). 


This problem has been resolved since 
long. Don't forget to activate Options > 

Simplify before Plotting or Options > Ap- 
proximate before Plotting in the Plot Win- 
dow. Otherwise you will receive an error 
message: 

2. Wir berechneten den Kreisinhalt. (We calculated the area of a circle): 



#1: radius(r, n, k) := 4 


" 2 

" k ■ r " 

2^ 

r - 




^ n ^ 

j 


#2: area(r, n, k) 


2-radius(r, n, k) ■ r 


n 

#3: appr_area(r, n) := I area(r, n, k) 

k=l 




Derive - USER - FORUM 


D-N-L#10 


#4: Precision := Approximate 

#5: Notation := Decimal 

#6: [2 ■ appr_area(l , 10), 2 ■ appr_area(l , 100), 2 ■ appr_area(l , 1000), 2 ■ appr_area(l , 10000)] 

#7: [2 . 90451S326 , 3.120417031, 3.139555466, 3.141391477] 

2 

#S: 2 ■ appr_area(r , 10000) = 3 . 141391477 ■ r -SIGN(r) 

FL00R(l/n) 


#9: ]im appr_area(l, n) = 2«]"im n ■ I 

n-wo n-*0+ k=l 


2 2 

- k -n ) 


2 FLQ0R(l/n) 2 2 

#10: ]im appr_area(r, n) = 2-r ■ 5IGN(r) ■ ] im m I 4a - k -n ) 

n-w» n-*0+ k=l 


In 1993 the limit was 0. 


R. Wadsack, Vienna, A 

There is a little trick to accelerate both the Cursor in text windows and the Cross in Plot- windows. 

You could include this command into any batch file. 

mode con rate = 32 delay = 1 

Marko Horbatsch, North York, CAN 

Dear DERIVE'rs 

particularly dear MATH teachers, 

It is really quite refreshing to see how active high school math teachers are in Europe by trying to do 
more interesting material in the classroom using DERIVE. Here in North America (I spent 10 years in 
Toronto, Canada by now after doing all my studies in Frankfurt up to Ph.D. in theoretical physics) 
matters aren't quite as bright as far as student preparation for university is concerned. While doing 
some work on Computer Assisted Instruction (CAI) for university kids, I thought that some of this 
might be useful for high school students as well. 

I am working on a set of course notes that will hopefully be published in not more than a year's time. 
Here is a little demonstration of how to use DERIVE to demonstrate functions of a complex variable. 
The idea that functions like sin and cos can be analytically continued in the complex plane is not evi- 
dent even to some university math students. The lack of intuition about this should be easily overcome 
by plotting sequences of z- and f(z)-points in separate screens (screens 2 and 3 respectively). Discrete 
point plotting was preferred over a continuous plot to show how each point z maps onto its f(z). For 
plotting purposes the following DERIVE objects were created. 

line 1: CV(z) generates an [x,y] vector from a complex number; 

line 3: ZV(r,th) generates a circle of radius r that will contain a parameter th(eta) to allow parametric 
plots in th(eta) once the radius r is specified. 

line 2: defines our function of interest that has poles at z = ± i. 

The three choices of circles in the z-plane z = 0.5, 1, 2 produce the very intriguing images f(z) and 
show the blow-up near the pole locations. 
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I am convinced that your clever high school students will love this! This subject area has received 
attention lately with the interest in maps for producing fractal sets. 

We have also developed some stand-alone CAI programs for linear algebra, complex variables and 
Newtonian mechanics. They will run on any EGA/VGA type PC with 256 kB of memory. You can 
obtain the whole package by sending a cheque over $ 20 (or Eurocheque for DM 30) to cover distribu- 
tion costs. Unlimited copying/use of the programs is granted. 

Specify the media type (5.25"(360/1.2M) or 3.5"(720/1.4M)). 

If you want more detailed into first, just write to me: 

Prof. Marko Horbatsch 
Physics & Astronomy 
YORK University, 4700 Keele St 
NORTH YORK, Ont. M3 J 1P3 
CANADA 



E. Sawada, Hawaii, USA 

I once attended a DERIVE session where David Stoutemyer was presiding. He drew - see the graph - 
using parametric equations. I used his idea and extended it. I used vectors to initiate vertical and hori- 
zontal translations. The following program illustrates this translation. Aloha. 

VECTOR (VECTOR ( [ [COS (x) * (COS (x) +1) +k, SIN (2*x) +m] , 

[2+1/ 2*COS (x/ 6-2 . 4) +k, 1/2*SIN (x/6-2 .4) +m] , [3/2+k, 1/4+m] , 

[SEC (SIN (x) ) +m] ] , k, -6, 6, 3) ,m, -3, 3, 3) 


DNL: T ry this nice plot! 
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An idea came to me when I was working with translation. If we can translate about some given point, 
why not translate about some given function. The following program will illustrate how we can move 
a function f(x) about another function g(x) using parametric equations. 

Moving a hypocycloid of four cups about a parabola: 



VECTOR ( [SIN (A) A 3+t, COS (A) A 3+3-t A 2] ,t, -2, 2, 0.25) 


Solving Physics Problems in the Classroom 
with DERIVE 2.0, Part 2 

Marko Horbatsch, York, Ontario 

Another attempt to make the differential equation - eigenvalue problem more understandable to 
physics undergraduates is through the solution of the differential equation as an initial problem with 
the eigenvalue as a trial parameter. In particular, it is emphasized how only by the choice of discrete 
energy eigenvalues one can obtain localizable eigenfunctions. 

The stationary Schrodinger equation with a quartic potential is solved for a selected eigenstate 
using the Runge-Kutta algorithm supplied in one of the support files. Any other symmetric potential 
could be substituted instead of the quartic one and higher eigenstates can be solved for as well with 
this technique. The aim of the demonstration is to show the following: if one solves the equation as an 
initial-value problem starting at x = 0 with appropriate conditions for a symmetric (or alternatively 
antisymmetric) solution and treats the eigenvalue as a trial parameter, one can satisfy the requirement 
of normalizability, dictated by physics, only for very special energy values (the eigenvalues of the 
Hamiltonian). Schrodinger's equation in one dimension (cf. eq. (3)) is converted to a pair of first 
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order equations (in units where h = m= 1): 

<P=X » X=2(V(x)-E) (p. (5) 

This makes the problem solvable by a standard technique for systems of first order differential equa- 
tions. The coding of the equations is shown below, the functions (pandx are denoted by y and yp. 


# 1 : 


# 2 : 


ODEApprox imation . mth — Copyright (c) 1990-2003 by Texas Instruments Incorporated 

cl + ( 1 im p) + 2-(c2 + c3) 

V-HJ_ + c3 


RK_AUX3(p, v, u_, cl, c2 , c3) 


6 


#3: 

#4: 

#3: 

# 6 : 


RK_AUX2(p , v, u , cl, 

c2) := RK_AUX3(p, v, u_, cl, 

c2, 

] im p) 

v^u_ + c2/2 

RK^UXl(p, v, u_, cl) 

:= RKJ\UX2(p, v, u_, cl, 

] im 

P) 


V-HJ_ 

_ + 

cl/2 

RK_AUX0(p, v, vO, n) 

:= ITERATES (u_ + RK_AUXl(p, v 

, u_ 

, lim p) , u_, vO, n) 

V-MJ_ 

RK(r, v, vO, h, n) := 

RKJUJX0(h.APPEND([l], r) , v, 

vO, 

n) 


4 

#7: V(x) := 0.5-x 

#S: RHS_Cy , yp, et) := [yp, 2- (VO) - et)-y] 

#9: (RK(RH5_(y, yp, 0.3), [x , y, yp], [0, 1, 0], 0.1, 50))u[l, 2] 

#10: (RK(RHS_(y, yp, 0.6), [x , y, yp], [0, 1, 0], 0.1, 50))u[l, 2] 


The Runge-Kutta routine RK(r,v,v0,h,n) contains in its argument list the derivative of the 
solution vector (r), the vector of the dependent variables itself (v), the initial condition (vO), the fixed 
step-size (h) and the number of steps (n). Note that v contains the independent variable as well. 



Fig. 3 Numerical solution of the stationary Schrodinger equation for the ground 
state of a quartic potential as an initial value problem with the energy as a 
trial value. Upper sequence of points: E tl = 0.5, lower sequence: E t2 = 0.6. 
The correct ground state eigenfunction is nodeless, decreases to zero as- 
ymptotically and has an eigenvalue of E tl <E 0 < E t2 . 
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Derive 's Runge-Kutta algorithm provides an example of a task that can be broken up into a few 
sequential steps and is thus ideally suited for this environment. It also shows that within a symbolic 
language one can handle automatically exceptions, where a function itself may not exist at a point 
needed for the calculation, but a limit can be obtained and is sufficient for the calculation (see lines #2 
- #6). The main reason for the use of the LIM ( f ( x ) , x , a ) function which calculates the limit of/(*) 
as x approaches a , is, however, the lack of a command for substitution and LIM in its first attempt just 
tries that. 

Note: Recent versions of DERIVE provide a SUBST-function, so it is not necessary to use 
the LIM-function. 

cl + SUBST (p , v, u_ + c3) + 2-Cc2 + c3) 

#1: RK_AUX3(p, v, u_, cl, c2, c 3) := 

6 


Another subject area that is being included in modern physics curricula is that of discrete maps. 
Figure 4 demonstrates that DERIVE 2.0 represents an ideal tool in order to study these. The relatively 
simple logistic (or quadratic) map is used as an example to show how radically the behavior of a map 
can change as an order parameter (in this example r ) is varied. There are applications of this map in 
many cases of physics (a model of turbulence, voltage-current characteristics in certain semiconduc- 
tors, etc.) and other sciences [5]. Derive 's way of allowing one to look both at the numbers and graphs 
and the fact that one has to look at the numbers and see how they were generated allows one to im- 
prove the understanding of this currently very fashionable subject. 

For r - values up until * 2.9 the map is not too exciting: after some number of iterations the 
(n+l) st iterate of 

x n+ \ = r x n (1 -x n ) (6) 

essentially has forgotten what its starting value x 0 was. (This fact itself is worth observing in Derive). 
Fig. 4 shows as a function of r the last 10 iterates of a sequence with typically 50 points. Up to r « 2.9 
these last 10 iterates practically fall on top of each other, such that they are displayed as a single point. 
For 2.9 < r < 3.4 one finds such sequences show bistability for large «, i.e. independent of the starting 
value Xo they will oscillate between two values, which, again, are only determined by the value of r. 
What happens in the vicinity of r\ * 2.9 is called a pitchfork bifurcation and it can be seen that this 
phenomenon repeats itself at r 2 * 3.44 at both parts of the fork and the sequence will oscillate between 
four values. This behavior of further bifurcation continues repeating itself, but is happening faster as r 
increases. A study of the r n - values at which bifurcations occur leads to a critical value r«, beyond 
which no regularity can be observed (chaos) and a law which relates all bifurcation points to a funda- 
mental constant (the Feigenbaum number 5): 


const. 



( 7 ) 


From the generated sequences of points one can obtain an estimate of the Liapunov exponent, 
which can be defined as the limit of the expression given in line 4 as s — > 0 and n — > oo. The conver- 
gence of the exponent with respect to both parameters could be studied. The sign of the Liapunov 
exponent can be used as a measure of the fact whether the behavior is chaotic or regular. 
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The details of the implementation shown below provide examples of three ways to generate the 
iteration with the functions FUN (line 2), FUN1 (line 3) and FUNO (line 10) depending on whether one 
wants pairs of iteration index and function value, just the last iterate or simply a sequence of the iter- 
ates. The function BIFURC in line 13 takes a sequence generated for given r and provides a vector of 
points to plot the last 10 elements at given r. To produce figure 4 one has to generate sequences for the 
various r- values, simplify line 13 and plot. 


#1: 

ltnap(r, x) 

:= r-x-(l - x) 


#2: 

FUN C r , xO, 

n) := ITERATES^ + 1, lmap(r, v )j, 

v, [0, x0], n) 

#3: 

FUNlCr, xO, 

n) := (ITERATE^ + 1, lmap(r, v )j 

, v, [0, xO], n)) 


FUNlCr, x 0 + e, n) - FUNlCr, xO, n) ^ 


£ J 

liapCr, xO, n, e) := 



n 


#5 

#6 


" 0 

0.1 

1 

0.09 

2 

0.0819 


1 DIM(w) 

xbarr(w) := ■ I w 

DIM(w) i_=l i_ 

FUN(1 , 0.1, 20) 


#8: FUN1C1, 0.1, 20) 

#9: 0. 03207763711 

#10: FUN0(r, x0, n) := ITERATES(lmap(r, v) , v, xO, n) 

#11: w := FUNOCl, 0.1, 20) 

#12: w := [0.1, 0.09, 0.0819, 0.07519239000, 0.06953849448, 0.06470289227, 0.06051642800, 0.05685418994, 0.05362179103, 

0.05074649455, 0.04817128784, 0.04585081487, 0.04374851764, 0.04183458485, 0.04008445236, 0.03847768904, 
0.03699715648, 0.03562836689, 0.03435898637, 0.03317844642, 0.03207763711] 


BIFURC (w, r, n) := VECTOR (T r, w 1, i, n - 9, n) 
#13: L -U 


■BIFURC (w, 1, 20) 


' 1 0.04817128784 ' 

1 0.04585081487 

1 0.04374851764 


#16: 

BIFURC_(r, xO, n) := VECT0R(j 

>, 

(FUNOCr, :< 0 , 

n, M 

#17: 

VECTOR CBIFURC_Ca, 0.1, 

20), 


O 
l n 


#18: 

VECTOR (BIFURC_(a, 0.1, 

20), 

a. 

1.05, 4.95, 

0.1) 

#19: 

VECTOR (BIFURC_(a, 0.1, 

20), 

a. 

1.025, 4.975 

, 0.1) 

#20: 

VECTOR (BIFURC_(a, 0.1, 

20), 

a. 

1.075, 4.975 

, o.l) 


I added expressions #16 through #20 to make the plotting procedure easier, Josef 
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II 


■ ■ i 





2.5 2.6 2.7 2.S 2.9 3 3.1 3.2 3.3 3.4 3.5 3.6 3.7 3.S 3.9 4 4.1 4.2 4.3 

Fig. 4 Iteration of the logistic map as a function of the order parameter r. The display shows bifur- 
cations in the accumulation points as r passes through certain critical values. At 3.57 
chaotic behavior (loss of regularity) can be observed. 


Another important topic in the upper year undergraduate curriculum is the solution of partial 
differential equations. Most equations discussed in the classroom context are separable and thus re- 
duced to ordinary equations that can be handled with Derive. In particular, Fourier representations of 
the solutions are commonly used and Derive 's capabilities to show the result are given in the last 
example. To prepare the students for situations where this technique does not work one has to study 
numerical techniques. Finite difference methods are easily applied to elliptic problems, such as Pois- 
son’s or Laplace's equation and the resulting systems of linear equations can be solved in Derive. For 
time evolution problems various combinations of finite difference representations of the temporal and 
spatial derivatives lead to either explicit or implicit methods [6]. Explicit methods result in recursion 
relations for many coupled parameters that represent the value of the unknown function at discrete 
points in space-time. 

The example shown in figure 5 deals with the numerical solution of the diffusion equation in 
one dimension: 

J- f(x, t) = c ~y f(x, t) . (8) 

ct ox 

This equation has many applications in the sciences, e.g., fort he description of heat conduction. Here 
it serves merely the purpose of demonstration as to what can occur in the discretization by an explicit 
finite difference algorithm. Replacement of the time derivative by a forward difference and of the spa- 
tial second order derivative by a central difference approximation [6] results in the recursion relation 

( fu =Axi„ tj)): 


ftj + 1 - fu + + fi-u ~ 2 fij ) • 


( 9 ) 
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Here the independent variables x and t have been made dimensionless by appropriate scaling 
(0 < x < 1, t > 0) and the new diffusion constant r has also absorbed a At/ Ax 2 factor coming from the 
finite difference operators. After imposition of the usual boundary condition f$ = 0 and choice 
of an initial condition such as, e.g .,/ (x, t = 0) = x(l - x) the recursion (9) can be solved. 


#1: 

n ;= 5 







r 

1 1 

#2: 

xvek := 

VECTOR 


i , i , 0 , n + 1 




< n 

+ i J 

#3: 

U0(x) : 

= - x* (x - 

i) 


uinit : 

;= VECTOR(U0(xvek ), i , n + 2) 

#4: 




i 


UPLUS( 

u , i , n. 

J : 

= 


If i 

< n_ 



#5: 

ui(i + 1) 



0 

UMINUS(u , i) ;= 
If i > 1 

#6: ur(i - 1) 

0 


RHS_(u , r) := VECTOR(IF(i <n+2Ai>l, u + r.(UPLUS(u, i, n + 2) + 
#7: i 

UMINUS(u, i) - 2 • u ), 0), i p n + 2) 
i 

#8: solution(r, m) := ITERATES ( RHS_(u , r) , u, uinit, m) 

plot(w, j) := VECTOR(|"xvek , w 1, i, n + 2) 

#9: L i j , i J 

y2 - yl 

#10: cutl i ne(xl, yl, x2 , y2) := (x - xl) + yl 

x2 - xl 


#11: cutl i ne(xl, yl, x2, y2) := 


y2 - yl 


■(x - xl) + yl 


# 12 : 


k x2 - xl 

cutl(a, b) := cutline(a , a , b , b ) 

12 12 


■X(xl, x, x2) 


#13: VECTOR I cutl ((plot(sol ution(0 . 5 , 10), j)) , (plot(sol ution(0 . 5 , 10), 

J=1 i 


j)) ), j, 1, 10, 2 

i + 1 


#14: 


f 6 

/ECTOR I cutl ((plot(sol ution(l, 
J=1 


(plot(sol ution(l, 5), j)) 



The first four statements serve to define a mesh and to set up the initial condition on this mesh. 
Lines 5 and 6 define functions that will return the right hand and left hand neighbors of a vector ele- 
ment provided that the index i is not at the boundary in which case zero is returned. The special hand- 
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ling of the endpoints is not necessary for the implementation of (9) as given on line 7, as the boundary 
condition is incorporated there directly. It is merely included to demonstrate a simple use of the I F- 
function, which is a new addition in DERIVE 2.0. 

RHS_ (u, r ) defines the updating for the vector of solution points, that is it defines equation 
(9). solution ( r , m) generates the iteration for given constant r and number of steps m. 
plot ( w, j ) is a function which generates a vector of points representing the solution at time-level j. 
In order to plot straight line segments between these points the function outline is defined. It makes 
use of Derive 's intrinsic CHI -function which returns one for xl < x < x2 and zero elsewhere. Shown 
subsequently is a statement which upon simplification leads to a vector of functions that consists of 
straight line segments connecting the solution points. A single plot call to the result provides the 
graphs shown in the windows of figure 5. It is shown that multiple loops can be generated in Derive, 
but it is also clear that such deep nesting of calls require carefulness of the programmer. 


The plot windows show such attempted solutions to (8) with 5 steps and r = 1 (right window) 
and using 10 steps and r = 0.5 (left window) for the same time domain. It turns out that this instability 
has nothing to do with round-off, as this calculation remains in the rational numbers, which Derive 
handles exactly. The discrete solution points have been connected by straight lines for convenience. 
Parabolic interpolation would have been more consistent with the second order differencing in space, 
but would have required more sophisticated programming for display purposes. The instability is seen 
to not appear immediately, but only after the solution begins to deviate from the originally parabolic 
form. The reason for this behavior is due to the fact that initially the spatial finite operator does not 
produce any discretization error. 



Fig. 5 Numerical solution of the dimensionless diffusion equation in one dimension by an explicit 
finite difference scheme. The calculation is done in exact precision rational number arith- 
metic and displays an instability of the algorithm for r > 0.5. Window left: r = 0.5 and 10 
time-steps taken, solutions for every second time-step are displayed in descending order; 
window right: same time-steps displayed, but calculated within 5 steps (r = 1.0). 


Marko Horbatsch 
Department of Physics and 

Applied & Computational Mathematical Science Programme 
York University, Toronto, Ontario, M3J 1P3, CANADA 
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Physics Problems on the TIs 

It is not too difficult to solve the system of differential equations on the TIs - Runge-Kutta- 
method is implemented. What one has to do is to rename the variables: the derivatives are 
yl and y2 and the independent variable x has to be changed to t. 

In the DE-mode choose the settings as shown below and enter the two equations together 
with the initial conditions: 


r-‘i=]z 

F£t Y F3 T FH Y FSt Y FfiT YF? FU 1 

f GRAPH FORMATS > 

r | 1 | 


Coordinates .... RECT + 
Grid?.".' 1 . T'i", ; : : ; off* 

fixes 0N-> 

Leading Cursor. OFF* 

Labels OFF* 

Solution Method RK-* 

Fields FLD0FF + 

, CEnter=SFHJE) CESC=CfiNCEL^ , 



USE * AND + TD OPEN CHOICES 



| Fi 7 | friT F£t Y F3 Y 
|Zoon|Edit| 

fh Yf5-*T fat Y ff Y 

Hflll|style|flxes...| 

1 

JLOTS 

Plot l: 
tG=G. 
-■yl '=y2 
yi 1 = 1 

et) ■ yl | et = .6 


^y2'=2-(.5 t 4 - 
yi2=0 
y3'=B 
yi3= 
y4 1 = 


w3' <t>= 

MAIN RAD EXACT DE 



As one cannot plot simultaneously two solutions (for et = 0.5 and 0.6) one has to save the 
graphs and superimpose them. The graph of function 0.5x 4 can be added via F6. 



Recommended websites for treating DEs with the Tl: 

www.seq.etsmtl.ca/ti/edqraphe.pdf (in French!) recommended by Michel Beaudin 
www.math.armstronq.edu/ti92/de/de.html recommended by Nils Hahnfeld 

A short program shows the bifurcation on the Tl-graph screen: 


|Fi7 | iriY fet i 
|Control 

fFjT'i 

1/0 

■FHt'I 

Gar 

If ind...Tl"lodeT 1 


bifurct 1 7 u ? st ? xG ? n!> 

Prgn 

Local x ? r, i : ClrDrau 
For r, 1 ? u ? st 
xO+x 

For i,l,n 

If i2in-9 Then 

If abs<x)<lQG0!PtOn r ? x 
Endlf 

r*x*( l-x)->x 
EndFor 

EndFor: EndPrgpil 

MAIN RAD APPROX FUNC 
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Original screen shots from DERIVE Versions 1.22 and 2.51 



corm®: 'ilPHira Center Delete Help hour Options Plot Scale Ticks Uirriem 
Zooil 

Enter Option 

Cross x-im ir& y:0 Stale StlB00 y:100 


Stale irlBSO 


Derive Z£h-plot 
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Independent Repeated Experiments 

Ales Kozubik, Bratislava 

This paper is dedicated to my daughter Susane which was born at the same night when I 
was typing it. 

In this short paper I would like to present one of possibilities of using DERIVE in teaching the Pro- 
bability Theory. The lectures on Probability Theory have been presented in University of Economics 
in Bratislava since the school year 1990/91. During the years 1990/91 and 1991/92 the programs Eu- 
reka and MathCAD were used and since 1992/93 the DERIVE program is used. 

Very important topic of this lectures are the repeated random experiments. In this article I deal with 
the independent repeated random experiments i.e. if the probability of the random events remains the 
same at each of repeated experiments 

At first I define the function pnom ( k, n) 

n! 

pnom(k, n) := 

#1: DIMENSION (k) 

n k ! 

i=l i 

10! 

pnom([3 , 3 , 4], 10), = [4200, 4200] 

31.31.4! . 

which gives for any natural number n and any vector k= [k\, k s \ the polynomial number 

ft ^ n\ 

k x ,k 2 ,...,k s ) k x \k 2 ! ... k s ! 

At first I deal with n independent repetitions of an experiment with two possible results (we can 
denote them as success and failure). In order to compute the probability that in these n repetitions the 
success will be set in /-times the Bernoulli formula 

fn\ 

p(\-p) n ~ l 

w 

is used, where p is the probability of the success in one experiment. This formula is realized in 
function ir ( i , n , p ) : 

i n - i 

#3: ir(i, n, p) := C0MB(n, i).p .(1 - p) 

The next question which we can answer is the question: What is the most likely count of success 
during n repetitions of the experiment? The answer is given by formula 

/ g ({n + 1) • p - 1; (pi + 1) • p ) 

realized in the function 
mlc(n, p) := 

If (n + l).p - 1 = FL00R((n + l).p - 1) 

[ (n + l).p - 1, (n + l).p] 

FL00R((n + l).p) 



#4: 
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In case the endpoints of the interval are integers, the function mlc (n, p) gives the vector containing 
both of them and in the opposite case it gives the integer contained in the interval. 

The function 

LNR(<x, p) := 

If LN(1 - oO/LN(l - p) = FL00R(LN(1 - <x)/LN(l - p)) 

#5: LN(1 - <x)/LN(l - p) 

FL00R(LN(1 - (x)/LN(l - p)) + 1 

answers the least number of repetitions for success to set in at least one with probability greater or 
equal to a. That is the least integer greater or equal than 

i - 1 

#6: af (i , p) := p.(l - p) 

gives the answer to the question: What is the probability the success will be set in at first in the i-th 
repetition? 



DIM(k) 



solves the problem of independent repeated events in the case of a possible results of the experiment. 
In this function n is the number of repetitions, k = [k u k s ] is a vector such that k\+...+k s = n and 
p = \pl, is a vector such that p\+ ... +p s = 1. 


Now, let Ai, ..., As are the possible results of the experiment with probabilities p i, ...,p s respectively. 
The probability that in n repetitions of that experiment the results Ai, ..., A^ will be set in k u ..., k s - 
times respectively is given by formula 


k\ , k 2 , • • • , k s 


'Pl X Pl 2 - Ps° 


which is realized in the function irs ( k, n, p) . 


Using the functions defined above we will solve the following elementary problems from probability 
theory: 


We will throw a die 21 times. 

a) What is the probability that number six will be thrown six times? 

b) What is the most likely count of throwing number six in these twenty replications? 

c) How many times do we have to throw the die to reach at least one number six with 
a probability greater than 0.9? 

d) What is the probability that the numbers 1,2, 3, 4, 5, 6 will be thrown 1,2, 3, 4, 5, 6 times 
respectively? 
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pi 9 


#8: Notation ;= Decimal 

Answer to question a) 

# 9 : 

Answer to question b) 

# 10 : 


i r 


6 , 21 , 

6 J 


= 0 . 07548933698 


mlc 


21 , 

6 y 


= 3 


Answer to question c) 
# 11 : 

Answer to question d) 

r 


#12 : i rs 


#13: 


lnr 


0.9, 

6 y 


= 13 


[1. 2. 3, 4, 5. 6], 21, 


1 1 1 1 1 1 

p I p P p 

6 6 6 6 6 6 

0.00009359689105 


Let my try solving the problems in the “traditional” way (Josef): 


c 


#14: 


#15 : SELECT 


BINOMIAL_DENSITY 


6 , 21 , 

6 ) 


= 0.07548933698 


v = MAX 


TABLE 




BINOMIAI DENSITY 


x, 21, 


, x, 1, 21, 1 


1 2 

kk k 


6 J 


f 

f 1 1 

y\ 

v, TABLE 

BINOMIAL_DENSITY 

X, 21, 

, x, 1, 21, 1 


k 

l 6 J 

yy 


U2 


#16: 


[[3, 0.2312786059]] 




r 

r 

" 5 " 

x T 

#17: 

CEILING 

NSOLUTIONS 

k 

1 - 

k 

, 6 . 

= 0.9, x 

yy 


= [13] 


1 

" i ' 

2 

" 1 " 

3 

' 1 > 

#18: 21 COMB(20 , 2). 

6 

, 6 , 

•COMB (18, 3). 

, 6 > 

•C0MB(15 , 4). 

, 6 > 


•COMB (11, 


5)< 


6 


•C0MB(6 , 6). 


1 ^6 


= 0.00009359689105 
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Bivariate Normal Distribution 

Dr. G.Marinell, Innsbruck, A 


Die Zufal Is variable A hat eine univariate Normalverteilung mit dem Erwartungswert // und 
der Varianz a , wenn sie folgende Dichtefunktion besitzt: 


-{x-nf /[2a 1 ) 


f N ( X = 

V2^cr 

Die Werte der Verteilungsfunktion dieser Normalverteilung 


W(X<a) = F N (a) = \"_J N {x \ju,<J 2 )dx 
konnen in DERIVE mit dem Befehl 


NORMAL (z, m, s) 

bestimmt werden, wobei p = m und a = s ist. (siehe auch BOHM J., (1992), Normal Distri- 
bution & Inverse, DNL#8). 


Die zweidimensionale Zufallsvariable (X, Y) hat eine (nichtsingulare) bivariate Normalvertei- 
lung, wenn sie folgende Dichtefunktion aufweist: 

/in I >(X 2 , P S ) = 


1 

• exp 

2na x a 2 -y/l - p 1 

r \ / 


1 

V A 

x-Mi 

1 2 ( l -p 2 ) [ 

°"i v 


\ / \ 2 1) 


-2p- 


X-jU, 


°i y 


y-ju 2 


\ a 2 J 


+ 


y-p 2 


V a 2 J 


m 


# 10 : 


1 

*e 

2 

2-tt. 7(1 - p ) 



2 2 

+ y )/ (2 - (1 - p )) 




Plot for p = 0.3 


1 • 

1 

7T . 

2 

2-7T.7(l - 0.3 ) 


- (X - 




e 
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BINORMAL calculates the distribution moment of the bivariate Normal Probability function. 
M = (a - pi )/<j 1 , v = (b - h2)/ct2, r = ctI 2/(ct1 * ct2) 


VECTOR 


4—4 

l,z,l, +1 

0.5 


■VECTOR 


1 -z /2 f 7r-z+u 

z , e , NORMAL | - NORMAL 


7 ( 2-70 


7(1- r) 


7r - z+v 


7(l-r) 


#3: BINORMAL (u,v,r):=- 


VECTOR 


4—4 
l,z,l, +1 




0.5 


■VECTOR 


1 -z /2 ( 7r-z+u 

z, e , NORMAL 


7 ( 2 - 77 ) 


7r-z+u ^ f Jr-z 

-NORMAL 

7(l-r) J l 7(1- 


7r - z+v 
r) 


1 -z / 2 
-•e 


L 7(2-77) 


z , -4 , 4 , 0 . 5 


1,3 


1 -z /2 
-■e 


s. 7(2*77) 


z, -4, 4,0. 5 


1,2 


BINO calculates the distribution function with the possibility to determine the accuracy 
by the interval [m,n] and the increment s. 

(In BINORMAL these values are: m = -4, n = 4, s = 0.5) 


VECTOR 1 , z , 1 


n-m ^ 
— * 1 


■VECTOR 


1 -z / 2 

z, e , NORMAL 


7(2-77) 


7r-z+u ^ C 7r-z+v- 

•NORMAL 


L 7(1- r) 


7(l-r)~ 


#4: BIN0(u, v, r,m, n, s):=- 


C n-m 

VECTORl 1 , z , 1 , +ll 


■VECTOR 


1 -z / 2 

z, e , NORMAL 


7(2-77) 


7r-z+u ^ C 7i”-z+v- 

-NORMAL 


7(1- r ) 


7(l-r)- 


r 2 ' 

1 -z /2 


L 7(2-77) 


, z,m, n, s 


1,3 


-} 


1 -z / 2 
•e 


L 7(2-77) 


, z,m, n, s 


1,2 


Example: Let X normal distributed with mean = 40 and stddev = 2.5 and 
Y normal distributed with mean = 50 and stddev = 5.5. 
Correlationcoefficent r = 0.29. 


What is the probability of p(X < 38 and Y < 53)? 

r 33-40 53-50 


#5 : BINORMAL 


, ,0.29 

2.5 5.5 j 


# 6 : 


0.1769063132 


f 33-40 53-50 ^ 

#7: BINO , ,0.29, -5, 5, 0.1 

l 2.5 5.5 ^ 


#3: 


0.176911 6233 
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I've received some letters and a lot of material to be used in classroom from Keith Eames, Senior Lec- 
turer, Kingsway, Camden College. Here are the main parts of two letters, one worksheet, a question 
and— I hope so - a suitable answer. 


Part of the 1st letter: 

... Though I am interested in the Use of DERIVE as a teaching tool I feel I do not have the 
time to produce an English DERIVE manual, but thanks for the offer of help. If you feel that 
the material I have produced is suitable then please feel free to use it as you think fit. A men- 
tion in any literature where it is used would be appreciated. 


I enclose a copy of a recent assignment given to students on curve fitting population statis- 
tics using DERIVE. I will send you copies of any further development work on DERIVE 


(see next two pages!) 

The second letter: 

Dear Josef, 

I wonder if anyone in the DERIVE User Group can help me out with a solution to a problem 
involving transformation of experimental data to a linear form? 

I enclose a copy of a recent question given to students on curve fitting of experimental data. I 
was surprised that DERIVE would not solve it, with the FIT command. 


The population of a colony of insects is measured at 
daily intervals: 

Number of days 0 1 2 3 4 5 

Population 40 50 80 110 150 220 


It is believed that p and t are related by an equation of the form p = a-b*. Estimate the values 
of a and b. 

When I attempted to fit an equation of this form, DERIVE failed to give any answer. So when 
this did not work, I thought I could try the transformation ln(p) = ln(a) + t ln(b) to the data in 
the 6x2 matrix. I found that I could not take natural logs of column 2 in the data matrix. 

A similar problem was encountered in trying to fit a curve to data which was in the form of 
T = a.l n , where the values of a and n had to be found. 

The Physicists and Chemists within the College wish to be able to perform transformations 
on experimental data, where the data is in matrix form. Is it possible to take, for instance, the 
log or square root of a single column of data values? They do not really want to take the log 
or square root of every single element in a column of a 50 x 2 matrix. 

(see page 26!) 
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FOUNDATION CM3 MODULE: CENSUS ASSIGNMENT 

When using Derive save your work and print out the graphs. Include the graphs in 
your write up. 


Population of USA from 100 to 1970 


Year, T 

Population (in millions) 

1900 

75 995 

1910 

91 972 

1920 

105 711 

1930 

123 203 

1940 

131 669 

1950 

150 697 

1960 

179 323 

1970 

203 212 


(a) Use Derive to plot the given data. Comment on your plot. 

Enter the data with T = year values and Y = population values. Since the T-values 
are so large they can be transformed using the formula: 

f 7-1900 
10 

This transformation results in the year 1980 being equivalent to t = 8 and the corre- 
sponding year values: 


f=0, 1,2, 3,4, 5, 6,7 

(b) Fit various polynomials to the data starting with a linear fit then a quadratic fit 
etc. up to a polynomial of degree 7. 

A polynomial of degree 7 would be given by: 

Y = OCq + CC jf + CCjt 2 + OC-^ + OC^t^ + QTjf 5 + OC^t^ + (Xjt 7 

where ao, a- i, «2 ai are constants. 

Plot each polynomial and comment on the closeness of each curve to the actual data 
points. 
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(c) In each case use your polynomial equation to predict the population of USA in 
1980. Find the true value of the population of USA in 1980 and comment on the ac- 
curacy of the predictions. 


Fitting the linear model 

A linear model will be of the form: 


Y — CC g + CC\t 

You will be required to find the two unknown coefficients a 0 and a- i. 

Derive should give you ao = 71 .4434 and ai = 17.5084, with the estimate for the 
population in 1 980 as 21 1 .51 0 million. 

To estimate the population value in 1980 substitute t = 8 in the equation of the line of 
best fit. Use the <Manage> <Substitute> command in Derive to carry out the substi- 
tution. (Now simply Substitute via the menu bar.) 


Fitting the quadratic model 

A quadratic model will be of the form: 

Y = a 0 + a l t + a 2 t 2 

You will be required to find the three unknown coefficients a 0 , ai and a 2 . 

Derive should give you ao = 79.0375, ai = 9.91363 and a 2 = 1 .08498, with the esti- 
mate for the population in 1980 as 227.7786 million. 

As a check the coefficients of the polynomial of degree 7 are: 




f 78.1631 " 

<X| 


0.58341 

a 2 


11.2415 



1.38854 

a 4 


-3.24255 

oc 5 


0.990631 

a 6 


-0.118051 



v 0.005001 , 


(There are slight differences to the coef- 
ficients obtained by using recent ver- 
sions of DERIVE. Compare with the re- 
sults on page 27. Josef) 
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#1: data := 


#2: 


#3: 


0 40 

1 50 

2 SO 

3 110 

4 150 

5 220 

DATUM Cd_3 := [d_ul, LN(d_u2)] 1 

0 3 . 6SS379454 

1 3.912023005 

2 4.332026634 

3 4.700430365 

4 5 . 010635294 

5 5.393627546 


DATLN(data) = 


#4: FIT_AUX(d_) := FIT([t, a + t-b], DATLN(d_)) 

FIT_AUX(d_) 


#5: EXP_FUNCCd_) := p = e 


0. 3463003373 -t 


#6: 

EXP_FUNCCdata) = (p = 33. 33231493 -e 


3 


C0MP(d_, u, t) :i VECTOR Cl" d_ 

d_ 

, SUBSTCu, t, d. 

#7: 

L k.l 

k , 

2 




0 

40 

33.332 



1 

50 

54 . 29263594 


0 . 3463 -t 

2 

SO 

76.79376310 

#3: 

COMPCdata, 33.332-e ) = 

3 

110 

103.6344525 



4 

150 

153.6671039 



. 5 

220 

217.3673109 


k, DIM(d_)) 


#9: 


#10: census := 


1900 

75 . 995 

1910 

91.972 

1920 

105.71 

1930 

123.2 

1940 

131.66 

1950 

150.69 

1960 

179.32 

1970 

203.21 


-10 0. 01351599339- t 

#11: EXP_FUNC(census) = (p = 5.546037374-10 -e ) 
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Quadratic Regression 


#21: quadcens(t) 

#22: quadcens(t) 

#23 : quadcens(S) 


FIT( 


[t, f-t + g-t + h], [[o,i, 


2, 3, 4 , 5, 6, 7], census j,|2] '3 


1. 0S42S5714-t + 9. 913309523 - t + 79.03366666 

227.7744235 



See first the exponential regression line together with the quadratic regression: 




DATA 








cl 

c2 

c3 

c4 

c5 


1 

1900. G 

0.0000 

paaTOia 




2 

1910.0 

1 . 0000 

91.972 




3 

1920.0 

2.0000 

1G5.71 




4 

1930.G 

3.0GG0 

123.2G 




5 

194G.G 

4.0GG0 

131.67 




6 

195G.G 

5.0GG0 

15G.7G 




7 

196G.G 

6.0GG0 

179.32 





rlc3=75 .995 


MAIN RAD ARRRDK 

FUNC 




jf LOT:' i 

Plot 4: 
Plot 3: 
Plot 2: 
-'Plot 1 :L 


y 1=78. 857250551381 ■( 1 . 1447240458596)* 


J 


^y2=5. 542281 2687963 e -10 ■( 1 . 0136081 18941> 


■^y3=. 010842857142857 ■ x 2 + "40.21 1026190** 
y4=B 


y4 <x>= 


MAIN RAD ARRRDK 

FUNC 



As we have 8 pairs of data we can try to find the polynomial of degree 7 (without transforming the 
year values) which leads to an interesting graph: 


UpH 

■ ^ x 5 + dx 4 + ex 3 + fx 2 + gx + h + cens(x) 

Done 

■solve(cens(1900) = 75.995 and cens(1910> 
ii g = 36424027238. 1 and h= "9.715056 e12 

■ soloeCcensCG^ = 75. 995 and censt 1} = 91 ■ 

a =.032571 and b= -.828506 and c = 8.2* 

■ 4 2 + 74.8G6424502191 ■ x + 75. 995 + cens7Cx^ 


... -lQlx -9715055724494 -8*v4<x> 

MAIN RAD AFFRDK FUNC E/30 




r^lz^Il^ceWgraphlM^hlD^^I:? 1 


the polynomial 

MAIN 

RAD AFFRDH FUNC 


Same graph but xres =10 


The polynomial with transformed years 
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Von der Binomial- zur Normalverteilung 

W. Propper, Numberg 

(This contribution "From the Binomial- to the Normaldistribution" is written in German. So 
I will try to give short summaries at some points where I find that it could be useful for non- 
german readers. Editor) 


1. Die Binomialverteilung 

Jacob Bernoulli (1655 - 1705) wird die Formel 

k n - k 

#1: B(n, p, k) := CQMB(n, k)-p -Cl - p) 

zugesprochen. Mit ihr lasst sich die Wahrscheinlichkeit fur genau k Treffer in einer Bemoullikette mit 
der Lange n und der Trefferwahrscheinlichkeit p berechnen. Fur einfache Falle kann B (n, p, k) „zu 
FuB“ berechnet, bzw. in Tafelwerken nachgeschlagen werden. Bei „krummen“ Werten versagen diese 
Hilfsmittel in der Schule, und der Wunsch nach einer einfachen Naherungsformel liegt nahe. 

Zur Veranschaulichung der Binomialverteilung dient das Histogramm. In ihm wird eine Wahrschein- 
lichkeitsverteilung x —> JV(x) in der Form von Balken uber dem Argument x derart dargestellt, dass 
die Balkenflache gleich (bzw. proportional) der Wahrscheinlichkeit W(x) ist. 

Wir konstruieren deshalb zunachst Balken der Breite d und der Hohe h, die symmetrisch zu x auf die 
Abszissenachse gesetzt werden. Die Funktion BALK (x, d, h) leistet dies. 

-2 i— | 

■■1 

— i 1 1 1— i— l 1 1 

-1 12 3 4 5 

#3: IJjBSfWlWSI ---1 

(Im Plot-Fenster muss die Options-> Display > Points auf Connected eingestellt sein, damit die Balken 
mit ausgezeichneten Randlinien erscheinen.) 

(The students should plot BALK(x=/c, d= 1, h=B(n,p,k)\ e.g. n = 20, p = 0.5, k e {..9, ...12,..}.) 
Die Bedeutung des Konstrukts BALK fur (Derive-ungeubte) Schuler kann an wenigen einfachen Bei- 
spielen gezeigt werden. 

Um zum Thema uberzuleiten, lassen wir nun an der Stelle k Balken der Hohe B ( n , p , k ) und der 
Breite 1 nacheinander plotten. Anschauliche Bilder ergeben sich beispielsweise mit n = 20,p = 0.5 
und ke {9 .. 12}. Dabei legt schon der erste Balken eine Anpassung des Koordinatensystems nahe. Mit 
einer geeigneten Skalierung und Positionierung des Centers erhalt man eine passende Umgebung. 
SchlieBlich schalten wir noch uber Options > Display > Color den automatischen Farbwechsel aus, um 
Farbgebungsmoglichkeiten (sofern vorhanden) auf wichtigere Falle anwenden zu konnen. 

Die Schuler vermuten schnell, wozu die o.g. Plot-Ubungen nutzlich sein konnen. Deshalb leuchtet die 
Syntax der folgenden Funktion HIST1 leicht ein. Sie liefert ein Histogramm der Binomialverteilung 
B (n,p) . 


#2: BALK(x, d, h) := 


x - O.S-d 0 
x - O.S-d h 
x + O.S-d h 
x + O.S-d 0 
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#4: HISTl(n, p) := VECTOR (BALK (k , 1, B(n, p, k)), k, 0, n) 


|HIST1(10, Q . 35)| 



Durch geschickte Variation von n und p lassen sich Charakteristika der Binomialverteilung veran- 
schaulichen: 

Bei festem p und veranderlichem n sieht man das „Auseinanderflie6en“, der Histogramme. 

Festes n und variables p zeigt die Symmetric im Falle p = 0.5 und die Verschiebung des Maxi- 
mums in den Fallen p < 0.5 bzw. p > 0.5. 

SchlieBlich konstruiert man Verteilungen mit konstantem Produkt n p. 

(Um ertragliche Rechenzeiten zu erhalten, stellen wir die Option-Precision auf Approximate und die 
Option-Notation auf Scientific. 3 signifikante Stellen sind jeweils genug. Bei der oben vorgeschlage- 
nen Einstellung des Plot-Fensters kann man mit n bis zu Werten von 60 gehen. In jeder Gruppe wird 
das Aussehen der Plots farblich durch Option > Display >Color-Plot variiert.) 

(Change n and p to show the influence of the parameters for the distribution, (p = const. & 
n = var., p = var. & n = const., n p = const.). The last step - n p = const - prepares for stan- 
dardisation - the maxima remain on the same position. We transform the coordinates so that 
the y-axis will become an axis of symmetry.) 


2. Die Standardisierung 


Im letzten Schritt (konstantes np) wird der Weg zur Standardisierung vorbereitet, da hier die Lage der 
Maxima erhalten bleibt. Wir verschieben deshalb das Koordinatensystem so, dass die Maxima aller 
Histogramme an derselben Stelle liegen. Die Wahl dieser Stelle an sich ist willkurlich, aber wenn die 
Ordinatenachse Symmetrieachse (fur p = 0.5) wird, ist auch Astheten gedient. Die Verschiebung er- 
folgt also um das Produkt n p. 


#6: |Xn, p) := n-p 

#7: HISTZCn, p) := VECTOR (BALK (k - p(n, p) , 1, B(n, p s k)), k, 0, n) 

#3: HIST2CL0, 0.35) 




■ 0,3 





■ 0.2 

■ 0.1 








-5 -4 -3 - 

2 - 

1 

■- 0.1 


l 3 4 5 


Die Histogramme HIST2 mussen nicht in so vielen Varianten dargestellt werden, wie HIST1, denn 
sie zeigen das angestrebte und schon voruberlegte Resultat. 
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Schwieriger ist es, die Charakteristika der neuen Verteilung herauszuarbeiten. Es handelt sich um eine 
Wahrscheinlichkeitsverteilung, die jedemx = k- n-p die Wahrscheinlichkeit JV(x) = B(n,p , k) zuord- 
net. Mit der SUM-Funktion lasst sich leicht der Ausdruck 

n 

#9: ERW(n, p) := I (k - n-p)-B(n, p, k) 

k=0 


#10: ERWC10, 0.35) = 0 


bilden und mit Simplify vereinfachen. (Achtung: Option-Precision auf Exact umstellen). Entsprechend 
kann dann auf die Varianz eingegangen werden: 

n 2 

#11: VAR_(n, p) := I (k - n-p) -B(n, p, k) 

k=Q 

#12: VAR_(10 s 0.35) = 2.275 


Als Ergebnis lasst sich festhalten: HIST2 ist im Erwartungswert auf 0 normiert, jedoch die Varianz ist 
noch abhangig von der Wahl von n und p. (Interpretiert man die Varianz als ein MaB dafur, wie stark 
sich die Werte einer ZufallsgroBe um den Erwartungswert scharen, so lasst sich das mit HIST2 gut 
veranschaulichen.) 

(The expected value is 0 for each (n, p), but the variances differ. Using o(n, p) we create the 
function NORM. Before plotting NORM its attributes should be discussed.) 

Nun ist einleuchtend, dass zur weiteren Vereinheitlichung die Varianz benotigt wird. Wir betrachten 
jedoch die Standardabweichung cr. 

#13: d(n s p) := -/(n-p-Cl - p)) 

Mit Hilfe von o(n,p) bilden wir nun die Funktion NORM: 



f 

' k - p(n, p) 

1 

'N 

#14: NORMCn, p) := VECTOR 

BALK 

, 

, cr(n, p)-B(n, p, k) 

, k, 0, n 



v, cr(n s p) 

cCn, p) j 

J 


Vor dem Plotten sollten ihre Eigenschaften diskutiert werden: 

Bei den vorherigen Beispielen war immer cr> 1, sodass nun Diagramme zu erwarten sind, deren Bal- 
ken sich naher um diey-Achse gruppieren (( k - \ ±)l <j<k- p), schmaler (l/cr< 1) und hoher sind 
(cr - B > B ). Gut einleuchtend ist die Erhohung der Balken um den Faktor a ; wenn gleichzeitig die 
Breite um diesen Faktor gestaucht wird. 

Die Plots von NORM mit verschiedenen Werten von n und p zeigen sofort, dass das Ziel der Verein- 
heitlichung erreicht ist. Die Histogramme der standardisierten Verteilung sind weitgehend unabhangig 
von n und p . Mit den nachfolgenden Funktionen ERWW ( n , p ) und VARI ( n , p ) lasst sich die Ver- 
einheitlichung numerisch ausdriicken (und ggf. der Wunsch nach einem Beweis wecken). 



#15: NORM (10 , 0.35) 
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n k - p(n, p) 

#16: ERWWCn, p) := I B(n, p, k) 

k=0 cr(n, p) 


#17: 


#13: 

#19: 


VARI(n, p) := 


n 

I 


k - pCn, 


p) '2 


-B(n, 


k=0 'v cr(n, p) J 
[ERWWC10, 0.35), ERWW(35, 0.47)] = [0, 
[VARIC10, 0.35), VARIC35 , 0.47)] = [1, 


P, k) 

0 ] 

1 ] 


(erww (n, p) andVARi(n,p) show standardisation numerically) 


3. Die GaulJkurve 

Das nachste Ziel ist nun, von der diskreten Wahrscheinlichkeitsverteilung zu einer in R definierten 
Naherungsfunktion zu gelangen. Um hier besser experimentieren zu konnen (mehrere Plots von NORM 
sind wirklich unubersichtlich!), fuhren wir die Funktion PKTE ( n , p ) ein. Da im Plot-Fenster die 
Option-State noch immer auf Connected steht, ergeben sich mit n>20 schon sehr ordentliche (ab- 
schnittsweise lineare) Naherungen. 


#21: 

PKTEC30, 

0.35) 

#22: 

PKTE (50, 

0.1) 

#23: 

PKTE (100, 

0.7) 



(The next goal is to find an approximation function. As some plots of norm are difficult to be 
read on one screen so we use PKTE -Options > Display > Points > Connected. Which func- 
tions have graphs similar to these plots - as GLOI , GAU1 , ....? 

Fur die Frage, welcher einfache Funktionsgraph diese Punkte verbindet, lassen sich Experimente 
durchfiihren, die auf Funktionen basieren, welche aus der Analysis bekannt sind. Sie miissen symmet- 
risch zury-Achse liegen, ein Maximum bei (0 / -0.4) haben und im Unendlichen verschwinden. (Im 
Plot-Fenster sollte die Zeichenoption auf Connected No umgestellt werden, um unnotige Linien zu 
vermeiden.) Von den Schulem werden schnell Funktionen der Art GLOI und GAU1 genannt. Weitere 
Beispiele werden sich daran orientieren. 

0.4 

GLQl(x) := 

#24: 2 

1 + x 


2 

#25: GAUl(x) := 0.4-EXPC- x ) 
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Diese Experimente fiihren schlieBlich zur GauBschen Fehlerkurve oder kurz Gaufikurve. 

2 ^ 


1 


#26: GAUS(x) := 


7(2-77) 


-EXP 


2 J 


(Selbstverstandlich lasst sich dieser Schritt, besonders wegen V(27i), nicht ohne Hilfe des Lehrers tun. 
Er lasst sich jedoch mit 


#27: 


2 

- x /2 

e dx = ,/2 - ,/tt 


sehr elegant begriinden. 


Ein Plot der Funktion GAUS lasst sich nun mit weiteren Darstellungen der Funktion PKTE verbinden 
(wobei weiterhin die Punkte diskret gezeichnet werden sollten). Hier lasst sich gut verdeutlichen, dass 
die Naherungsfunktion gute Werte liefert, falls o> 3 ist. 


Fur eine Naherungsfunktion der Binomialverteilung muss der bisherige Weg zumckverfolgt werden: 
Das Argument der Funktion GAUS wird durch (x - ju)/c7 ersetzt und der Wert von GAUS wird durch 
den Faktor o{n,p) reduziert. Als Ergebnis erhalten wir 


#23: 1 

GAUSSfn n — -F V P 

f 

r x - p (n, p) ^ 

2 ^ 

tj(n, p) > 

'JiThU OiJ ^1 1 ^ |w* p ■ — LAP 

cj(n, p)-7(2-7T) 

2 j 


(At last we have to consider the transformation x — > (x - p)/a and the normation - divide by 
a (n, p), so we obtain #28.) 


Eine Wahrscheinlichkeitsverteilung, deren Dichte durch die Funktion GAUSS bestimmt ist, heiBt nor- 
mal verteilt. Vergleiche der Plots von GAUSS mit HIST1 (bzw. der nachfolgenden Funktion PKTB) 
zeigen, dass das Ziel erreicht ist. Entsprechend konnen numerisch die Werte von B (n, p, k) und 
GAU S S ( n , p , k ) verglichen werden. 


#29: pktb(n, p) := VECTGR([k, B(n, p, k)], k, 0, n) 


#30: HIST1CL00, 0.35) 

#31: pktbO-00, 0.35) 

#32: GAUSS (100, 0.35) 



(Finally we compare the plots of GAUSS with HIST1 and PKTB and the numerical values of 
BIN(n,p,k) with GAUSS(n,p,k).) 
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3. Praktische Anwendung 


Zum Schluss ware noch die Frage anzuschneiden, wie sich nun die in der Stochastik haufig auftreten- 
den Problemstellungen "Berechne die Wahrscheinlichkeit fiir eine Trefferanzahl im Intervall \k\ \ k2 ]’’ 
oder auch ’’Berechne P(X < £)’’ bearbeiten lassen. 

Uber die Binomialverteilung liefert die Funktion WBIN eine Losung (wobei dies in der schulischen 
Praxis ohne Computerunterstiitzung mit der ebenfalls tabellierten kumulativen Binomialverteilung 
angegangen wird; kritischen Riickfragen kann durch Vergleich von WBIN(n,p,0,k) mit 
O ( n , p , k ) begegnet werden. 

k2 

#33: WBIN(n, p, kl, k2) := I B(n, p, k) 

k=kl 

Aus der Erinnerung an die Herleitung des Riemannschen Integrals kommt schnell die Vermutung, dass 
bei Verwendung der Funktion GAUSS ein Integral der Form 

k2 

#34: WGAl(n, p, kl, k2) := f GAUSS(n, p, x) dx 

kl 

das Problem losen wird. Die Berechnung von WBIN und WGA1 mit einigen konkreten Wertebelegun- 
gen zeigt aber, dass die erste Vermutung falsch war. Insbesondere ist immer WGA1 < WBIN. Ein Ver- 
gleich mit den letzten Plots zeigt die Ursache auf: 

#35: WBIN (100, 0.35, 32, 40) = 0.641S672045 
#36: WGA1(100, 0.35, 32, 40) = 0.5330622500 

Bei WBIN sind die Randstreifen voll beriicksichtigt, wahrend sie bei WGA1 nur etwa zur Halfte mit 
eingehen. Deshalb ist die in WGAU angebrachte Korrektur erforderlich. 

k2 + 0.5 

#37: WGAU(n, p, kl, k2) := J GAUSS(n, p, x) dx 

kl - 0.5 

#33: WGAU (100 , 0.35, 32, 40) = 0.6440325550 

Die Aufgabe ’’Berechne P (X < £)” fiihrt demnach auf die Funktion 
_k + 0.5 

#39: 4(n, p, k) := J GAUSS(n, p, x) dx 

—CO 

#40: WBIN (100, 0.35, 0, 42) = 0.940569S166 

#41: $(100, 0.35, 42) = 0.9420742512 

Hiermit ware diese Darstellung zu einem gewissen Abschluss gebracht. 

(We want to show that there is a good approximation for p(k1 < X < k2) when we use the 
integral but use the limits kl-0.5 and k2+0.5). 
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Arrows and Labels for the Axes 

Josef Bohm, Wiirmla 

(Later versions of DERIVE allow annotations in all plot windows. So it is no more necessary to design 
one’s own x and y to label the axes. Arrows might be still useful.) 


Letters X and Y 






m - 0. 05 -sx 

n - 0. 05 -sy 





0.05-sx + m 

0. 05 -sy + n 

#1: 

LX (m, 

n, 

sx, sy) := 

m 

n 





m - 0.05-sx 

0. 05 -sy + n 





_ 0.05-sx + m 

n - 0. 05 -sy . 





m - 0.05-sx 

n - 0. 05 -sy 





0.05-sx + m 

0. 05 -sy + n 

#2: 

LY(m, 

n, 

sx, sy) := 

m 

n 





_ m - 0.05-sx 

0. 05 -sy + n . 

The arrow heads up / down / right / left 






m - 0.04-sx 

n - O.OS-sy 

#3: 

AUOn, 

n, 

sx, sy) := 

m 

n 





_ 0.04-sx + m 

n - O.OS-sy . 





m - 0.04-sx 

n + O.OS-sy 

#4: 

AD(m, 

n, 

sx, sy) := 

m 

n 





_ 0.04-sx + m 

n + O.OS-sy . 





m - O.OS-sx 

n + 0.04-sy 

#5: 

AR(m, 

n, 

sx, sy) := 

m 

n 





_ m - O.OS-sx 

n - 0.04-sy . 





m + O.OS-sx 

n + 0.04-sy 

#6: 

ALOn, 

n, 

sx, sy) := 

m 

n 





_ m + O.OS-sx 

n - 0.04-sy _ 


sx and sy are the values for the scales 
in x- and y-direction - or multiples 
for bigger arrows and letters. 

(m, n) is the location of the center of 
the letters and the peaks of the ar- 
rows. 


2 

x 


#25: 

y = 4 

2 


#26: 

AU(0, 5, 2, 4) 


#27: 

AR(4, 0, 2, 4) 


#2S : 

LX(4.2, 0.3, 1, 

2) 

#29: 

LYC0.3, 5.2, 1, 

2) 



